Comment on "Critical current density from magnetization 
hysteresis data using the critical-state model" 

P. Chaddah 
Low Temperature Physics Laboratory, 
Centre for Advanced Technology, 
Indore 452013, India 
(February 6, 2008) 

Abstract 

A recent paper [Physical Review B 64 014508 (2001)] claims to present 

an exact method to extract the critical current density Jc from the M-H 

hysteresis curve. We show that this claim appears unjustified. 
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The hysteresis in the isothermal M-H curves of superconductors was first related to the 
critical current density 3c by Bean's critical state mo defl. Fietz and Webbi showed that 
the errors involved in extracting the field-dependent Jc(B) from these M-H curves are of 
order (d 2 Jc*/dB 2 ). The applicability of this method of extracting Jc(B) has been studied 
in great detail!! after the discovery of high-T^ superconductors. The recent claim! of an 
'exact method' of obtaining the same information from the same data, consequently, needs 
to be examined. 

The 'exact method' of reference [5] relates, in its equations (7) and (8), the derivatives of 
the forward and reverse envelope M-H curves to Jc(B a ), and to the difference between the 
applied field B a and the field B c + (or B c _ ) at the centre of the sample during the increasing 
(or decreasing) field cycle. It then presents a method (see figures 2 and 3 of reference[5]) to 
relate the central fields B c + (or B c _) to B a via the measured M-H curves. All these results 
are obtained for a sample in the shape of an infinite slab parallel to the applied field. 

Analytical solutions of Bean's critical state model, for various Jc(B), are known for 
infinite cylinders (of various cross-sections) in parallel field geometry, when the sample de- 
magnetization factor N=0 .!H! It is known that while field profiles B(x) depend on Jc(B) 
and not on the shape of the cylinder's cross- sect ion, the value of the magnetization M is 
dictated by the cylinder's cross-section .!! While equations (5) and (6) of reference[5] are 
valid for all N=0 cylinders, the derivation leading to equations (7) and (8) breaks down. 
This is simply because the volume element at x depends on the shape of the cylinder's 
cross-section, resulting in different x-dependent weight-functions within the integral sign. 
Equations(ll) of reference[3] give explicit analytic expressions for the M-H curves of circu- 
lar cylinders when Jc(B) follows the exponential model Jc(B) = Jc(0)exp(-B/H ). It is also 
known! for the exponential model that B^ + = H log[exp(B a /H ) - (H*/H )] and that Bf_ 
= H log[exp(B a /H ) + (H*/H )]. It is straightforward to confirm that equations (7) and 
(8) of reference[5] are not satisfied. 

The argument leading to figure 3(b) of reference [5] is based on the relation (depicted in 
figure 2 thereof) that the total flux contained in the sample is the same in the field-increasing 
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case when the applied field is B a , as in the field-decreasing case when the applied field is 
E>c+ (corresponding to that B a ). It follows from the results in reference[3] that this is not 
true for circular cylinders for Jc(B) following the exponential or the Kim- Anderson models. 
It follows from reference [1] that this is not true for a circular cylinder even when Jc(B) is 
independent of B. It is straightforward to see that this result also breaks down for cylinders 
with other cross-sections!. 

We assert that the method of reference [5] is applicable only for an infinite slab in parallel 
field, and is not valid for any other sample shape. The conventional method^Hi for extracting 
Jc(B) from isothermal M-H hysteresis curves is valid for various shapes, and involves errors 
whose upper bounds are known. 

Finally, reference [5] claims to show that the shape of the hysteresis bubble identifying the 
peak-effect in Jc(B) is asymmetric about the peak field Bp. This is shown using the input 
that Jc(B) decays as the magnitude of (B-Bp) rises. While this result has been shown earlier 
explicitly with reference to the peak effect (see figure 2b of reference [7]), the calculation is 
identical to the special case of Bp=0. This asymmetry in the 'central peak position' has 
been studied extensively in recent years!. 
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